Running coupling: Does the coupling between dark energy and dark matter change 

sign during the cosmological evolution? 



Yun-Hc Lii and Xin Zhangi'^,* 

^Department of Physics, College of Sciences, Northeastern University, Shenyang 110819, China 
^Center for High Energy Physics, Peking University, Beijing 100080, China 

In this paper we put forward a running coupling scenario for describing the interaction between 
dark energy and dark matter. The dark sector interaction in our scenario is free of the assumption 
that the interaction term Q is proportional to the Hubble expansion rate and the energy densities 
of dark sectors. We only use a time-variable coupling 6(a) (with a the scale factor of the universe) 
to characterize the interaction Q. We propose a parametrization form for the running coupling 
b[a) — boa + be(l — a) in which the early-time coupling is given by a constant be, while today the 
coupling is given by another constant, bo- For investigating the feature of the running coupling, we 
employ three dark energy models, namely, the cosmological constant model {w = —1), the constant 
w model (w = wq), and the time-dependent w model {w{a) — vuq + wi{l — a)). We constrain the 
models with the current observational data, including the type la supernova, the baryon acoustic 
oscillation, the cosmic microwave background, the Hubble expansion rate, and the X-ray gas mass 
fraction data. The fitting results indicate that a time-varying vacuum scenario is favored, in which 
the coupling b{z) crosses the noninteracting line (b = 0) during the cosmological evolution and the 
sign changes from negative to positive. The crossing of the noninteracting line happens at around 
2 — 0.2 — 0.3, and the crossing behavior is favored at about la confidence level. Our work implies 
that we should pay more attention to the time-varying vacuum model and seriously consider the 
phenomenological construction of a sign-changeable or oscillatory interaction between dark sectors. 

PACS numbers: 95.36.-l-x, 98.80.Es, 98.80.-k 



I. INTRODUCTION 

Dark energy (DE) and dark matter (DM) are the dom- 
inant components in the current universe, according to 
the recent astronomical observations [1-3]. They to- 
gether account for about 96% of the critical energy den- 
sity of the universe. However, ironically, we have known 
httle about the natures of DE and DM. Although DM is 
"dark" (nonluminous) , its gravitational property is nor- 
mal, i.e., it gravitationally behaves like the usual baryon 
matter and thus can form structures in the universe. The 
property of DE is much more exotic in that it is gravita- 
tionally repulsive and so responsible for the accelerated 
expansion of the universe [4]. 

Since we are ignorant of the natures of both DE and 
DM, we cannot ignore such an important possibility that 
there is some direct, non-gravitational interaction be- 
tween DE and DM. Intriguingly, such a possible interac- 
tion between DE and DM plays a crucial role in helping 
solve (or, at least alleviate) several important theoretical 
problems of DE. For example, it can be used to under- 
stand the cosmic coincidence problem [5], to avoid the 
cosmic doomsday brought by phantom [6], and to solve 
the cosmic age problem caused by old quasar [7] as well. 
Therefore, it is very meaningful to seriously study the 
interaction between dark sectors. Owing to the lack of 
the knowledge of micro-origin of the interaction, one has 
to first phenomenologically propose an interacting DE 
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model and then test its theoretical and observational con- 
sequences. So far, lots of phenomenological interacting 
DE models have been studied [8, 9]. 

When considering the interaction between dark sec- 
tors, the continuity equations for energy densities of DE 
and DM are of the form: 

Pde + iH{l+Wde)pde = ~Q, (1) 
Pdm + SHpdm = Q, (2) 

where pde and pdm are the energy densities of DE and 
DM, respectively, H = a/a is the Hubble parameter, a 
is the scale factor of the Friedmann-Robertson- Walker 
(FRW) universe, Wde = Pde/ Pde is the equation of state 
(EOS) parameter of DE, a dot denotes the derivative 
with respect to cosmic time t, and Q denotes the phe- 
nomenological interaction term. Several forms for Q 
have been put forward and have been fitted with ob- 
servations [10]. Of course, all of these models are phe- 
nomenological. Most of them are constructed specifically 
for mathematical simplicity — for example, models in 
which Q (X Hp, where p denotes the energy density of 
the dark sectors, and usually it has three choices, namely, 
P = Pdm, P ^ Pde, and p = Pde + Pdm- In addition, there 
arc also some models [11] in which the assumption of that 
Q is proportional to the Hubble parameter is abandoned 
and thus Q cc p. Such models are designed by consult- 
ing the simple models of reheating, of dark matter decay 
into radiation, and of curvaton decay — i.e., where the 
interaction has the form of a decay of one species into 
another, with constant decay rate. However, it should 
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be stressed that these models are severely dependent on 
the man-made choice of the special interaction forms. 
In other words, the predictions and observational conse- 
quences are model-dependent. In particular, the above- 
mentioned phenomenological models exclude the impor- 
tant possibility that the interaction changes sign during 
the cosmological evolution. It is of interest to point out 
that a sign-changeable or oscillatory form of interaction 
is possible, according to the current observations. 

Recently, Cai and Su [12], without choosing a special 
phenomenological form of interaction, proposed a novel 
scheme in which the whole redshift range is divided into 
a few bins and the interaction term S{z) (note that here 
Q = 3H6) is set to be a constant in each bin, and by fit- 
ting the observational data they found that S{z) is likely 
to cross the noninteracting (6 = 0) line. This study is 
fairly enlightening and inspires us to open our mind to 
seriously consider the possibility of that the interaction 
between dark sectors changes sign during the cosmolog- 
ical evolution. If such an observation is conclusive, it 
is suggested that more general phenomenological forms 
of interaction term should be put forth. However, the 
work of Cai and Su [12] seems not sufficient to prove that 
the interaction changes its sign, since their conclusion is 
drawn based upon the behavior of the best-fitted S{z), 
and the errors of the fitting results weaken the conclusion 
to a great degree due to the fact that the observational 
data currently available cannot determine more than two 
parameters in a piecewise parametrization approach. 

In this paper, our aim is to verify whether the inter- 
action indeed changes its sign (i.e., crosses the nonin- 
teracting line) during the evolution, by using a different 
method. We propose a parametrization form for the in- 
teraction term Q. In our work, we further abandon the 
assumption that Q is proportional to the Hubble expan- 
sion rate H. So, the interaction term Q is only charac- 
terized by the coupling b: 

Q{a) = Sb{a)HoPo, (3) 

where the dimensionless coupling 6(a) is variable with 
the cosmological evolution. Note that here the occur- 
rence of the present-day Hubble parameter Hq and the 
present-day density of dark sectors po ~ pdeo + Pdmo is 
only for a dimensional reason. By the way, in the whole 
work the subscript "0" always indicates the present-day 
value of the corresponding quantity. From the form of 
Eq. (3), we see that the evolution of the interaction term 
Q is totally described by the running of the coupling con- 
stant, 6(a), so our scenario can be called the "running 
coupling". Furthermore, we assume that the coupling b 
is described by a constant bo at the late times, and de- 
termined by another constant be at the early times; and 
the whole evolution of 6(a) is totally characterized by the 
two parameters, 6o and be- For continuously connecting 
the early-time and late-time behaviors, we put forward 
the following two-parameter form for the coupling 6(a): 

6(a) = 6oa + 6e(l-a). (4) 



Though the interaction term in our work depends on a 
particular parametrization form, the parameters can be 
tightly constrained by the current observational data, 
overcoming the disadvantage of the piecewise fitting 
method. The reconstructed evolution of 6(a) will indicate 
whether the coupling between the dark sectors crosses the 
noninteracting line. 

In this paper, we will investigate our coupling 
parametrization with the latest observational data. For 
the EOS of DE, w, we consider the following three 
cases: (1) the cosmological constant (vacuum energy), 
w = —1; (2) the constant EOS, w = wq; (3) the time- 
variable EOS, namely, the Chevallier-Polarski-Linder 
(CPL) parametrization, w{a) ^ wq + wi{l — a) [13]. We 
will fit the three interacting DE models with the data 
from the Union2 type la supcrnovae (SNIa), the baryon 
acoustic oscillation (BAO), the cosmic microwave back- 
ground (CMB), the Hubble expansion rate, and the X- 
ray gas mass fraction. We obtain the best-fitted param- 
eters and likelihoods by using the Monte Carlo Markov 
chain (MCMC) method. We will show that the interac- 
tion Q between DE and DM indeed changes sign around 
z = 0.2 — 0.3 during the cosmological evolution, at about 
la confidence level (CL). 

II. METHODOLOGY 

We consider interacting DE models in a spatially fiat 
FRW universe. The Friedmann equation reads 

MIpiH^ = Pr + Pb + Pde + Pdm, (5) 

where pr, Pb, Pde and pdm are the energy densities of 
radiation, baryon, DE and DM, respectively, and Mpi is 
the reduced Planck mass. It is convenient to introduce 
the fractional energy densities 51^ = pi/3AIpiH^ , with 
i = r, 6, de, and dm. Obviously. 

fir + ^b + ^de + ^dm = 1- (6) 

Substituting Eqs. (3) and (4) into Eqs. (1) and (2), 
and defining the functions fde = Pde/pdeo and fdm = 
Pdm/Pdmo, we obtain 

^^+3{l+w)fde{x) = (^1 + i.^ [5oe- + 6,(1 - e^)] , 

(7) 

+ ^fdrnix) = t|^(1 + /o) [&0e" + foe(l ~ e^] , 

ax ^yx) 

(8) 

where a; = Ina, /o = Pdeo/Pdmt^ = ^deo/^dmo = 
(1 - Vtro - VLbo - ^dmo)/^dmO, and E{x) = H{x)/Ho = 

[flroe^"^^ + flboe^^^ + ^dmofdrn{x) -f (1 — flrQ — ^bO ~ 

^dmo).fde{x)Y^^ ■ Therefore, given the values of the pa- 
rameters flro, ^bo, ^dmo and w, Eqs. (7) and (8) can be 
numerically solved with the initial conditions fdeiO) = 1 
and /(jm(0) = 1. With the resulting functions fde{x) 



3 



and fdm{x), we finally obtain the function E{x), the di- 
mensionless Hubble expansion rate. As aforementioned, 
in this work we employ three DE models, namely, the 
cosmological constant model (ACDM) with w = —1, 
the constant EOS model (XCDM) with w ^ wq, and 
the time-variable EOS model (CPL) with w{x) ~ wq + 
wi (1 — e^). 

To fit the three interacting DE models with observa- 
tions, we use the data from the Union2 SNIa (557 data), 
the BAO from SDSS DR7, the CMB from 7-year WMAP, 
the Hubble expansion rate (15 data), and the X-ray gas 
mass fraction (42 data). The best-fitted parameters are 
obtained by minimizing the sum 
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We obtain the constraints by using a MCMC method. 

Supernovae. — We use the data points of the 557 
Union2 SNIa compiled in Ref. [14]. The theoretical dis- 
tance modulus is defined as 



= 51ogiQ DL{zi) + ^iQ, 



(10) 



where z = l/a— lis the redshift, = 42.38 — 5 log^Q h 
with h the Hubble constant Hq in units of 100 km/s/Mpc, 
and the Hubble-free luminosity distance 



Dl{z) = {l + z) 



dz' 



E{z'-e) 



(11) 



where denotes the model parameters. Correspondingly, 
the function for the 557 Union2 SNIa data is given by 



[y-obsjZi) - y.th{zi)Y 



(12) 



where a is the corresponding Icr error of distance modu- 
lus for each supernova. The parameter [Iq is a nuisance 
parameter and one can expand Eq. (12) as 



x|^(0)-^(0)- 2^0^(0) 



(13) 



where A(6»), B{B) and C are defined in Ref. [15]. Evi- 
dently, Eq. (13) has a minimum for /ig = BjC at 



c 



(14) 



Since xlw, 



XsN, 



, instead minimizing Xsn 



will minimize Xsn which is independent of the nuisance 
parameter /ip. 

Baryon acoustic oscillations. — We use the BAO data 
from SDSS DR7 [16]. The distance ratio (d^) at z = 0.2 
and z = 0.35 are 



do.: 



^0.35 



(15) 



where rg{zd) is the comoving sound horizon at the baryon 
drag epoch [17], and 



Dv{z) 



dz' 



H{z') J H{z) 



1/3 



(16) 



encodes the visual distortion of a spherical object due to 
the non Euclidianity of a FRW spacetime. The inverse 
covariance matrix of BAO is 



(C 



BAo) ^ 



30124 -17227 
-17227 86977 



(17) 



The x^ function of the BAO data is constructed as: 



xlAO = {df-df-){C^\oMd 



'th 



jobs\ 



(18) 



where di ^ ((io.2j 1^0.35) is a vector, and the BAO data we 
use are do. 2 — 0.1905 and do. 35 — 0.1097. 

Cosmic microwave background. — We employ the 
"WMAP distance priors" given by the 7-year WMAP ob- 
servations [18]. This includes the "acoustic scale" /a, the 
"shift parameter" i?, and the redshift of the decoupling 
epoch of photons z» . The acoustic scale I a describes the 
distance ratio D a{z<,) / r s{z<,) ^ defined as 



U = (1 + 2* 



^t:Da{z* 
rs{z*) 



(19) 



where a factor of (1 4- z*) arises because Da{z^,) is the 
proper angular diameter distance, whereas ^^(z*) is the 
comoving sound horizon at z,. The fitting formula of 
rs{z) is given by 



1/(1 + 2) 



da 



a2i?(a)v/l + (3f^bo/41^7o)a 



(20) 



In this paper, we fix fi^o ~ 2.469 x 10 ^, 17bo ~ 
0.02246/1^2, given by the 7-year WMAP observa- 
tions [18], and rj^o = ^^70(1 + 0.22717Veff) with N^s the 
effective number of neutrino species (in this paper we 
take its standard value, 3.04 [18]). We use the fitting 
function of z* proposed by Hu and Sugiyama [19] 

z, = l{m[l + 0.m\2A{nMh^)-^'^^\l + qi{Vl,nahy% 

(21) 

where Vtmo = ^ba + ^dmo and 



2\-0.238 



9i 



0.0783(r26o/i2) 

l + 39.5(f7fco/l2)0-763' 



.92 



0.560 



l + 21.1(f7fco/l2)l-81' 

(22) 

The shift parameter R is responsible for the distance ratio 
Da{z*)IH-^{z,), given by [20] 

i?(z,) = ^n,noH^{l + z,)Da{z.). (23) 

Following [18], we use the prescription for using the 
WMAP distance priors. Thus, the x^ function for the 
CMB data is 

XcMB ~ {^T ~ ^f"'){^c\lB)i'A^*j' ~ ^j^*)' (24) 

where Xi = {Ia,R,z,:) is a vector, and (C'cMB)«i 
inverse covariance matrix. The 7-ycar WMAP observa- 
tions [18] give the maximum likelihood values: lAiz*) = 
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302.09, R{z^) = 1.725, and = 1091.3. The inverse 
covariance matrix is also given in [18]: 



i^CMB ) 



2.305 29.698 -1.333 
29.698 6825.27 -113.180 
-1.333 -113.180 3.414 



(25) 



Hubble expansion rate. — For the Hubble parameter 
H(z), there are 15 observational data available, where 
twelve of them are from Ref. [21]. In addition, in 
Rcf. [22], the authors obtain the additional three data: 
Hiz = 0.24) = 79.69 ± 2.32, H{z = 0.34) = 83.8 ± 2.96, 
and H{z = 0.43) = 86.45 ± 3.27 (in units of km s'^ 
Mpc~^). The function for the observational Hubble 
data is 
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i=l 



[Hthje; z,) - Hobsiz,)f 



(26) 



X-ray gas mass fraction. — For the fgas data, we use 
the Chandra measurements in Ref. [23]. In the frame- 
work of the ACDM reference cosmology, the X-ray gas 
mass fraction is presented as [23, 24] 



fgas{z) 



n 



60 



KAjb{z) 

1 + s{z) V ilbo + ndrnoF{z) 



n 1.5 



Da{z) 



(27) 

where the effect of the interaction between dark sectors 
has been considered, resulting in an additional function 
F{z) = fdni{- ln(l + z))/{l + zf. The parameters K, 7, 
b{z) and s{z) model the abundance of gas in the clusters. 
We set these parameters to their respective best-fit values 
of Ref. [23]. A is the angular correction factor, which is 
caused by the change in angle for the current test model 
^2500 in comparison with that of the reference cosmology 

nACDM . 
"2500 



V ^2500 / 



H{z)DAiz) 
[H{z)DAiz)]^^^''^ 



(28) 



here, the index r] is the slope of the /gas(?'/''250o) data 
within the radius ^2500, with the best- fit average value 
rj = 0.214 ± 0.022 [23]. And the proper (not comoving) 
angular diameter distance is given by 



Da{z) 



1 



dz' 



{l + z)Jo Hiz') 



(29) 



The function for the fgas data from the 42 galaxy 
clusters reads 



Xx- 



rayW 



^ {[fgasjO; Zi)]th 
i=l 



[fgasiZi)]obs) 



<r^{z^) 



(30) 



It should be pointed out that the fgas data are rather 
crucial for the fitting, since they can be used to break the 
degeneracy between the parameters from the interaction 
and the EOS of CPL model. So, the inclusion of the fgas 
data in our fitting is indispensable. 



III. RESULTS 

Now we fit our interacting models with the obser- 
vations. For the interacting ACDM, XCDM and CPL 
models, the parameters are 6 = {fidmo, ^e, h}, 
{17dmO, Wo, 60, &e, h) and {^drao, Wo, wi, bo, be, /i}, re- 
spectively. We use the MCMC method and finally we ob- 
tain the best- fit parameters and the corresponding Xmin- 
The best-fit, Icr and 2cr values of the parameters with 
Xmin of the three interacting models are all presented in 
Table I. 
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FIG. 1: The probability contours at la and 2a confidence lev- 
els in the parameter planes for the interacting ACDM model. 



Figure 1 shows the likelihood contours for the inter- 
acting ACDM model. For this model, we have ildmo — 
0.2262, bo = 0.0793, b^ = -0.3274 and h = 0.7120, 
with Xmin = 595.968. We plot the likelihood contours 
for the interacting XCDM model in Fig. 2. For the in- 
teracting XCDM model, the fitting results are ildmo = 
0.2267, Wo = -0.9844, 60 = 0.0787, be = -0.3216 and 
h =^ 0.7094, with Xm^n = 595.815. For the interacting 
CPL model, we obtain the fitting results: fldmo = 0.2271, 
Wo = -0.9768, wi = -0.0455 60 = 0.0799, be = -0.3290 
and h = 0.7097, with Xm^n = 595.808. The likelihood 
contours for this case arc shown in Fig. 3. 

From the fitting results, we first notice that the EOS of 
DE is near —1 in both the constant w and time-dependent 
w scenarios. In the constant w case, we get the best-fitted 
value Wo = —0.9844, rather close to —1, albeit slightly 
tends to a quintessence {w > —1). For the CPL case, the 
best-fitted values are wo = —0.9768 and wi = —0.0455; 
we also notice that wo is fairly near —1 and wi is close 
to 0. So. in the interacting DE model with a running 
coupling, even if the EOS of DE is allowed to vary, the 
data still favor a slightly evolving EOS with the value 
approaching —1. To sec this clearly, we plot the recon- 
structed w{z) for the interacting CPL model in Fig. 4 
where the best fit and the Icr uncertainties are shown. 
According the fitting results, we realize that a time- 
varying vacuum scenario is favored by the data. The 
time-varying vacuum model has been discussed exten- 
sively [25, 26]. Our work indicates that the scenario of 
a time- varying vacuum with a running coupling deserves 
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TABLE I: The fitting results of tlie parameters witli best-fit values as well as la and 2a errors in the three interacting DE 
models. 



model parameters 
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XCDM 


CPL 


^dmO 
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Wl 

bo 
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FIG. 2: The probability contours at la and 2a confidence levels in the parameter planes for the interacting XCDM model. 



more further investigations. 

For the running coupling 6(z), wc obtain similar results 
in all the three scenarios. From Figs. 1-3 wc sec that 
the parameters bo and be are in strong anti-correlation. 
The best-fitted values for feg and be are: bo ~ 0.08 and 
be ~ —0.3, implying that the coupling b{z) crosses the 
nonintcracting line 6 = during the cosmological evo- 
lution, and the sign changes from negative to positive. 
Such a feature for the interaction is favored by the data at 
about 1(7 level. To see the crossing phenomenon clearly, 
we reconstruct the evolution of the coupling b{z) by using 
a Fisher Matrix technique, shown in Fig. 5. From this fig- 
ure we read out that the crossing happens at z = 0.2—0.3. 
Our results tell such a story: at early times when DM 
dominates the universe, the energy transfer direction is 
from DM to DE, and at late times when DE becomes 
dominant, the decay direction reverses, from DE to DM. 
The above phenomenon is favored at Icr CL. Note that 
once the reconstruction of b{z) is performed based on a 
Monte Carlo method, the fluctuations will become larger, 
but the above distinctive feature still stands by at about 
la CL. 

In the work of Cai and Su [12], a redshift-binned 
parametrization method is used, but the limitation is 
that it is hard to go beyond two parameters for tight 
constraints. Also, in Rcf. [12], as the data give only 
weak constraint for z > 1.8, the additional assump- 
tion Q = for 1.8 < z < 1090 is made. In our 



method, the parametrization Q{a) = 3b{a)HoPo with 
b{a) = boa -I- foe (1 — a) has only two parameters, and the 
early-time interaction can also be described. Our fitting 
results support the conclusion of Cai and Su [12] that 
the interaction Q{z) crosses the nonintcracting line, in 
a distinct way. The limitations of our parametrization 
are: (1) whether or not there is some oscillation feature 
in the interaction cannot be read out, and (2) the future 
evolution cannot be described owing to the fact that b{z) 
will diverge as z approaches —1, so the predictive power 
of this parametrization is lost (see [27] for a similar case 
for w{z)). We will go beyond these limitations in the 
future work. In addition, now that the analysis of the 
current observational data provides a hint that the inter- 
action between dark sectors might change sign during the 
cosmological evolution, more general phenomenological 
forms for the interaction describing the sign-changeable 
or oscillatory feature should be seriously considered. 



IV. CONCLUSION 

Since the knowledge about the micro-origin of dark 
sector interaction is absent, one has no way to know the 
form of the interaction term Q from a microscopic theory. 
A popular way to investigate the interacting dark energy 
is to assume a specific phenomenological form for Q; for 
instance, Q on Hp or Q cc p, where p denotes the energy 
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FIG. 3: The probability contours at la and 2a confidence levels in the parameter planes for the interacting CPL model. 
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FIG. 4: The reconstructed evolutionary history for w{z) of 
the interacting CPL model. 



density of dark sectors. However, the choice of the phe- 
nomenological forms is rather arbitrary. In the face of 
this situation, let us recall the method used to probe the 
dynamical evolution of the EOS of DE, w{z), with the 
observational data. Owing to our ignorance of DE, for 
probing the dynamical evolution of w, one has to param- 
eterize w empirically, usually using two parameters, e.g., 
the most widely used CPL form, w{a) = wq + 'Wi{l — a) . 
Inspired by this method, we try to parameterize the evo- 
lution of Q in a similar way, and then use the observa- 
tional data to probe the dynamical evolution of Q. Such 
a way may provide a guidance for finding a reasonable 
phenomenological form for Q. 

In this paper, we have put forward a running cou- 
pling scenario for describing the interaction between DE 
and DM. In this scenario, the dark sector interaction 
has the form Q{a) = 36(a)iJoPo, where b{a) is the cou- 
pling which is variable during the cosmological evolu- 
tion. We have proposed a parametrization form for the 
running coupling: b(a) = boa + &e(l — a). So, at the 
early times the coupling is given by a constant be, while 
today the coupling is described by another constant, 
6o- We have constrained the parameters bg and b^ with 



the observational data currently available, including the 
Union2 SNIa, BAO (from SDSS DR7), CMB (from 7- 
year WMAP), H{z), and X-ray gas mass fraction data. 
It should be stressed that the /gas data are very crucial in 
our analysis since they are helpful in breaking the degen- 
eracy between the parameters. And, in our analysis, we 
employed three DE model, namely, the ACDM, XCDM 
and CPL models. 

Our fitting results show that the EOS of DE, w, is 
close to —1 both in the constant w (XCDM) and time- 
dependent w (CPL) cases. Thus, a time- varying vacuum 
scenario is favored by the data, according to this analy- 
sis. In addition, for the running coupling b{z), the results 
are also similar in all the three DE models. The param- 
eters 6o and be are in strong anti-correlation, and the 
best-fitted values are: bo « 0.08 and be ~ —0.3. This im- 
plies that the coupling b(z) crosses the noninteracting line 
{b = 0) and the sign changes from 6 < to & > 0. The re- 
construction of b{z) indicates that the crossing of the non- 
interacting line happens at around z = 0.2 — 0.3, and the 
crossing behavior is favored at about la CL. Therefore, 
our work tells a story about the interacting DE model: 
DE is a time-varying vacuum; the coupling between DE 
and DM runs with the expansion of the universe; at early 
times when DM dominates the universe, DM decays to 
DE, while at late times when DE becomes dominant, DE 
begins to decay to DM. If the above scenario is true, 
then we should pay more attention to the time-varying 
vacuum model, and seriously consider how to construct 
a sign-changeable interaction between dark sectors phe- 
nomenologically. For the time-varying vacuum model, 
previous work neglects the perturbations of DE and only 
treats the decaying vacuum as a background, however, if 
DE interacts with DM, the perturbations of DE should 
also be taken into account even if u; = —1. This deserves 
further investigations. 

Finally, we discuss the limitations of our paramctriza- 
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tion for the running coupling. Our paramctrization can- 
not describe the possible oscillation in the dark sector 
interaction. Moreover, such a paramctrization form can- 
not predict the future evolution of the interaction, since 
the coupling b(z) will encounter divergency in the far fu- 
ture. We will go beyond these limitations in our future 
work. 
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